Absence of stable collinear configurations in Ni(OOl) ultrathin films: canted domain 

structure as ground state 
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Brillouin light scattering (BLS) measurements were performed for (17-120)A-thick 
Cu/Ni/Cu/Si(001) films. A monotonic dependence of the frequency of the uniform mode on an 
in-plane magnetic field H was observed both on increasing and on decreasing H in the range (2- 
14) kOe, suggesting the absence of a metastable collinear ground state. Further investigation by 
magneto-optical vector magnetometry (MOKE-VM) in an unconventional canted-field geometry 
provided evidence for a domain structure where the magnetization is canted with respect to the per- 
pendicular to the film. Spin wave calculations confirm the absence of stable collinear configurations. 
PACS numbers: 75.70.Ak, 75.70.Kw, 75.30.Gw, 75.30. Ds, 78.35. +c, 78.20.Ls 



I. INTRODUCTION 



The determination of the magnetic phase of thin films 
can require a subtle analysis owing to the presence of 
many different competing anisotropics. Furthermore, 
metastable phenomena can be present, so that the mag- 
netic history of the sample must be considered. For ex- 
ample, when the magnetization is perpendicular to the 
film plane, a monodomain state is suggested from rema- 
nence of polar Kerr loops, El even thougJji a striped do- 
main configuration is the ground state.H However, it is 
worthwhile to note that, in general, all the investigated 
films present either perpendicular (0 = 0) or in-plane 
[6 = 7r/2) magnetization (6 is the angle between the 
magnetization and z, the normal to the film plane); a 
phase with canted (0 < 8 < 7r/2) magnetization has been 
usually observed close to the reorientation transition be- 
tween perpendicular and in-plane magnetization in a very 
nacraw range of film thickness, t, and/or temperature, 
Tan This phase is usually related to a domain structure 
as experimentally observed by microscopic techniques 
Canting of the magnetization has been observed in a class 
of rather different systems as well: multilayered struc- 
tures of Fe/Tb, for example B 

It is generally agreed that epitaxial Cu/Ni/Cu/Si(001) 
films present a collinear phase magnetized perpendic- 
ularly to the film plane ip._^i_svide range of Ni thick- 
ness, 17 A< t <120 A.M£3 Assuming the perpen- 
dicular state to be stable, the frequency f^o of the uni- 
form mode versus H, an external magnetic field applied 
within the film plane,._is predicted to have a minimum 
at a critical field i/ c .lid In a previous paper ,E3 we re- 
ported Brillouin light scattering (BLS) experiments in 
Cu/Ni/Cu/Si(001) films which, in the thickness range 17 



A< t < 80 A, showed a monotonic decrease of the fre- 
quency of the uniform mode on decreasing the in-plane 
field down to 2 kOe, with a clear change in slope at field 
values H* w (5-8) kOe, depending on the film thick- 
ness. (A similar behavior was.-abserved also in other, 
systems: Co/Au ultrathin filmsjij Co/Pd multilayers,!! 4 ] 
and Co/Au multilayers.ta) We attributed such a discrep- 
ancy between our data and theoretical predictions! 1 !! to 
the breakdown of the collinear ground-state state and to 
the onset of a domain-coafiguration.E£l As a matter of 
fact, spin- wave theorylail 8 ] predicts such an instability 
at a field value H> = H c + AH C (see Section V for the 
definitions of H c and AH C ) and at a finite wavevector 
fey , whose inverse sets_the initial size of the domains. 
In our previous paperO strong evidence for the pres- 
ence of domains was provided by magneto-optical vector- 
magnetometry (MOKE-VM) data. A direct visualization 
of up and down perpendicular domains in these speck 
mens for zero applied field was presented by Hug et aL;Ej 
unfortunately, our present attempts to visualize - using 
magnetic force microscopy - the domain structure while 
applying an in-plane magnetic field were unsuccessful due 
to the very weak magnetic signal from the sample and to 
the interference with the strong external field. 

If the system is saturated in the z (perpendicular) di- 
rection and the in-plane field H is increased from to 
H c , a collinear ground state with a canting angle 6 in- 
creasing from to 7r/2 is expected, and the spin-wave 
frequency is predicted to decrease until, at a field value 
Hf Ri H c — AH c /2, a "lower" instability (quite similar 
to the "upper" one at H>) is expectedO The fact-that, 
on decreasing H, the Brillouin light scattering datalij did 
not show any indication for the low field phase might be 
interpreted as due to a metastability phenomenon. 
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The aim of the present paper is to definitively inves- 
tigate the possible occurrence of the above mentioned 
metastability and to achieve a detailed understanding of 
the ground state and of both the static and dynamic mag- 
netization properties. New BLS data taken by increasing 
the in-plane field H are presented and discussed. More- 
over, MOKE-VM data in a new canted-field geometry are 
presented in order to analyze the collinear ground state. 



II. BRILLOUIN LIGHT SCATTERING DATA 



Our new measurements of BLS spectra are reported in 
Fig. 1 for the Cu/Ni/Cu/Si(001) sample with Ni thick- 
ness t = 60 A for several values of the external magnetic 
field H applied within the film plane (similar results are 
obtained for other thicknesses in the range 17A-80A). 
The details of the film growth and structural charac- 
terization have been published elsewhere, E2I but a brief 
description is necessary here. The films were grown in 
a molecular beam epitaxy (MBE) chamber by e-beam 
evaporation. The base pressure was less than 2.0 10~ 10 
Torr and rose in the low 10 -8 Torr range during deposi- 
tion of both Nickel and Copper layers. The growth rates 
for the Copper and Nickel layers were 1.0 A/s. The cop- 
per buffer layers for all the samples were annealed in situ 
to about 150°C for eight minutes to improve the flat- 
ness of the buffer layer surface. Unlike the case of Cu 
on Si(001), the interface betwen Ni(001) and Cu(001) re- 
mains sharp at_coom temperature. It has been shown 
by Chen et al.cil that the Ni on Cu is thermally sta- 
ble against interdiffusion up to 200° C. Reflection high 
energy electron diffraction (RHEED) and synchrotron x- 
ray diffraction show high quality epitaxial Nickel films. 
BLS measurements were carried out in the backscatter- 
ing geometry at the GHOST laboratory in PerugiaE3 us- 
ing a Sandercock-type tandem Fabcy-Perot interferome- 
ter in a (3+3)-pass configurations About 150 mW of 
P-polarized light, from an Ar + -ion laser operated in sin- 
gle longitudinal mode on the 5145-A line, was focused 
onto the sample surface, at an incidence angle of 45°, 
by a camera objective of numerical aperture 2 and focal 
length 50 mm. The sample was placed between the poles 
of an electromagnet used to produce a d.c. magnetic field, 
with a maximum intensity of 15.0 kOe, applied parallel 
to the film surface and perpendicular to the plane of in- 
cidence of light. Since light scattered by magnons has its 
plane of polarization rotated through 90°, an analyzer 
was used to remove unwanted back-reflections from the 
objective lens and light scattered by acoustic phonons. 
The spectra, recorded at room temperature, were stored 
in 512 channels, each with a gate length of 1 ms. A few 
thousand interferometer scans were necessary to record 
spectra with a good signal-to-noise ratio. 



III. MODEL 

From the measured spectra in Fig. 1, one obtains the 
field dependence of the frequency of the Damon-Eschbach 
mode reported in Fig. 2. A detailed theoretical analysis 
of these data is performed in terms of the following spin 
Hamiltonian 



H = 



i,S i i 

+ W E^) 4 + \ K l E^f) 4 + (^) 4 ] 
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, (Si ■ rjj)(Sj ■ Tjj) 



(1) 



where J > is the nearest neighbour ferromagnetic ex- 
change constant; H is an external field; K2 > is a uni- 
axial anisotropy favouring the z direction, normalJx) the 
film; K^r > and k\ are quartic anisotropics fcnH-J the 
former favours the film plane, xy; the latter, when the 
spins are in plane, favours them to lie along the in-plane 
diagonal [110]; finally, the last term in Eq. (1) is the 
magnetostatic dipole-dipole interaction. The presence of 
a quartic in-plane anisotropy is required since the BLS 
data, performed for field applied in plane along the easy 
[110] and the hard [100] direction respectively, show a 
small difference in frequency (see the inset in Fig. 2). For 
sufficiently high fields (H > H>), the spins are aligned 
in plane along the field direction and the frequency of the 
isQ 



uniform mode 



" I [ i j .oj = 7\/ (H - H K ott - -H k k ) {H + H K \\ 



(2a) 



tto|[ioo] = lJ{H - H Kf - H K ,){H - „) (2b) 



for H along the easy and hard in-plane axis, respectively 
(7 is the gyromagnetic factor). Here above and in the 
following we put H K »tc — Hk 2 — Hdi P , where g^BHK 2 = 
IK2S is the uniaxial anisotropy field and Hdi P = AttMqCi 
the dipolar field (M is the bulk magnetization and 
ci = //\/2 is a thickness dependent coefficientOEJ) ; 
g(j, B H \\ = 21<ls 3 and g/j. B H K ± = 2K^S 3 arc the 
quartic anisotropy fields. The former is estimated to be 
H K \\ « 1 kOe from the frequency difference in the BLS 

data, (Au w (2-3) GHz for H > 10 kOe), while a fit of 
flo(H) at high fields along [110], using Eq. (pa]), yields 
an estimate for the effective anisotropy field H K «tt w 2.5 
kOe. Finally, we define the exchange field giJ,BH ex = 2JS 
which - though not appearing in Eqs. ( |2a| , |2b| ), since they 
are written for zero wavevector - will be useful in the 
following. 

The results of our new BLS measurements in Fig. 2 
clearly show that almost the same curve f2o versus H is 
obtained on increasing as well as on deceasing the in- 
plane field H in the range (2-14) kOe.EZI We conclude 
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that the collinear perpendicular state is a very weakly 
metastable state or that it is plainly unstable, i.e. the 
system does not present a perpendicularly magnetized 
phase. Even more important, BLS data exclude a ground 
state with domains magnetized perpendicularly to the 
film plane as well. In fact, for H = and negligible 
quartic anisotropics, the energy gap of the spin- wave ex- 
citations with respect to a ground state with stciped do- 
mains perpendicularly magnetized is estimated!^! to be 



1 -L, stripe 



1 {H K2 - H dip N zz ) 



(3) 



where N zz < 1 is the demagnetization factor. Thus, in 
a perpendicular striped domain structure, the gap would 
be higher than in the case of spin waves excited from a 
uniform perpendicular ground state (i.e. with N zz = 1) 



\ uniform 



(4) 



Since, using the value of H K ^a deduced before, the latter 
frequency is estimated to be at least 7 GHz, the hypothe- 
sis of a perpendicular striped domain structure in the low 
field region (H < 2 kOe) appears to be inconsistent with 
the high field (H > 10 kOe) results. This problem can 
be solved assuming the system to have a quartic perpen- 
dicular anisotropy K± > 0, i.e. competing with K^* and 
leading to a canted collinear configuration for H = 0. Ne- 
glecting the small in-plane quartic anisotropy, the zero- 
field canting angle 8o is given by cos 2 0q — H K df / H K ± 



(see Eq. (A4) in the Appendix). In this way, the problem 
of the inconsistency between low and high field BLS data 
is removed because, for H — > 0, the spectrum of the ex- 
citations with respect to such a canted collinear ground 
state has a Goldstone mode (Qq = 0) as a consequence of 
the rotational invariance around the z axis. In the pres- 
ence of a small quartic in-plane anisotropy, the gap at 
zero field is slightly greater than zero, but small enough 
to make our argument still valid. 



IV. MOKE-VM DATA 

In order to ifiyeal the possible occurrence of a canted 
magnetization!^ we use the magneto-optic Kerr vector- 
magnetometry (MOKE-VM) techniqueEJ in a novel 
canted-field geometry. The field is applied within the 
plane formed by z = [001], the normal to the film, and 
x = [110], the in-plane easy axis; the field direction 
formed a canting angle Oh with z. In Fig. 3a,b we re- 
port the MOKE-VM data for the polar (M z ) and the 
longitudinal (M x ) component of the magnetization, re- 
spectively. In the present experimental configuration we 
did not detect any transversal component (M y ) in the 
whole investigated field range. In Fig. 3c, we show the 
field dependence of the calculated magnetization modu- 
lus M = \J Af 2 + M 2 , normalized to the saturation value 
M sat . 

The different curves refer to different orientations of 
the applied magnetic field, with Oh ranging between 0° 



and 20°. The perpendicular component M z shows a re- 
manence in zero field, the greater the lower is Oh- Also 
the parallel component M x is found to be nonzero in zero 
applied field; moreover, when the applied magnetic field 
is sufficiently high (\H\ > 4 kOe), M x is found to de- 
crease with increasing H for Oh < 11° and to increase 
for Oh = 20°. Therefore, these observations call for a 
canted ground state configuration in zero field with a 
canting angle 0q comprised between 11° and 20°. For 
the canted-field g eom etry of the MOKE-VM experiment 
one has (see Eq. (A4) in the Appendix) 



cos 6> = (H K e 



\h k1 )/(H k , + ±H k i). (5) 



Hence, putting 0q ~ 15°, we obtain H K ± w 2.7 kOe. 
Concerning the magnetization modulus, we observe that 
it is always decreasing as the field is decreased, mainly 
owing to the the decrease of the M z component. This is 
a strong evidence for a slightly canted ground state with 
a domain structure for the perpendicular component of 
the magnetization: M = (M x , 0, ±M Z ). 



V. DISCUSSION 



Now, after having estimated all the Hamiltonian pa- 
rameters which determine the frequency of the uniform 



mode (i.e. H K v 



H 



2.5 kOe, H 



1 kOe and 



K ± ~ 2.7 kOe) we note that, using these values, the 
BLS data are well reproduced (see Fig. 2) only in a very 
high field region (H > 10 kOcL. Lel order to discuss the 
possible arising of instabilities,E£rE£l we need to general- 
ize the calculation of the spin-wave frequency gap to a 
finite wavevector, f2o(k||). For doing that, we specialize 
to the in-plane easy axis direction [110]. 

For long wavelength, the spin-wave acoustic mode 
(which is identified with the Damon-Eschbach mode in 
the Brillouin spectrum) of an ultrathin film with N 
planes takes the approximate form 



n (k||)« 7 xMi(k||) MH) 



(0) 



where the full expressions of A\ (k|| ) and A2 (k\\ ) are given 
in the Appendix. This expression is approximated in the 
sense that it is the generaliza tioji n£.the spin-wave disper- 
sion relation of the monolayeirni 16 l'Ll±o the acoustic mode 
of an ultrathin film with N planes.EJ Since Ai(kii) is al- 
ways strictly positive, the only reason for a field-induced 
instability at finite wavevector can be the vanishing of 
A2 (k|| ) = ao — ai (fey a) + a2(k\\a) 2 (where a is the square 
lattice constant and the oti 's are known functions of the 
physical parameters, see Eq. (A2b) in the Appendix). 



The minimum of ^(kii) is obtained for k 



011/(20:2) and we obtain ^(fen™-*) — ao — a^/^o^): if 
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this quantity is negative, the spin wave frequency be- 
comes pure imaginary, signaling the instability of the uni- 
formly magnetized canted state. The instability appears 
when As (fey ) — 0: this condition defines the values 

of the upper (H>) and lower (Hf) critical fields. In 
the high field regime we find = H c + AH C , with 



H 



and a 



1 TJ 

i n K\ 



and AH r 



(H dip N) 2 / (16 f 2 H^), 
the uniform ground 
s 1/fc 



(m) 
I 



(H dlp N)/(AfH ex ) 

state breaks into domains of size 

The BLS data show the existence of two relevant fields: 
at H rs 9 kOc experimental data start deviating from 
the theoretical spin wave frequency and at H* ~ 6 kOe 
there is a clear change in their slope. The higher field is 
interpreted to signal the arising of a domain structure, 
which makes the experimental fio deviate from the value 
predicted for a collinear in-plane phase: it is therefore 
identified as the field H> . Instead, the changement in 
the slope at a field H* « 6 kOe should signal that a well 
defined domain structure is now established, whose size 
may strongly interfere with the light used in the BLS 
experiment. 

We can sum up our interpretation as follows. At high 
fields we have an in-plane collinear phase and we start 
decreasing H. At H = H> ss 9 kOe, the collinear 
phase is unstable against the appearance of a domain 
structure: spins acquire a very small z component that 
changes sign from a domain to a neighbouring domain. 
This instability corresponds to the vanishing of (k|| ) 

for a value which is just the inverse of the domain 
size. With further decreasing H, the domain structure 
becomes more and more defined and the domain size in- 
creases. 

On the other side, the "lower" instability occurs for 
w H c — AH c /2. In our present interpretation, the 
quantity AH C is taken as a free parameteio and esti- 
mated, from the high field results, to be ss 6 kOe. Con- 
sequently, since H c ps 3 kOe, the "lower" instability is ex- 
pected to occur already for H sa 0. We observe that such 
a conclusion about the instability of a collinear canted 
ground state is strongly supported by our MOKE-VM 
results in canted field geometry: even a relatively high 
field H sa 4 kOe applied nearly along the easy magne- 
tization axis (Oh = H° or 20°) is unable to give the 
saturation value of the magnetization modulus, see Fig. 
3c. 

In conclusion, new BLS measurements, performed 
upon increasing the in-plane field H, allowed to def- 
initely exclude a mctastable collinear state for epitax- 
ial Cu/Ni/Cu/Si(001) films with thickness in the range 
(17-80)A. This feature is confirmed by spin-wave calcu- 
lations. The MOKE-VM data in canted-field geometry 
suggested the presence of a domain structure where the 
magnetization is slightly canted (0 ss 15°) with respect to 
the film normal. We are confident that our results can be 
useful to explain spin wave anomalies previousljf-ebserved 
in other systems, such as Co/Au ultrathin filmstU, as well 



as Co/Pd multilayersEi and Co/Au multilayers. 



APPENDIX A: SPIN- WAVE DISPERSION 

The spin-wave dispersion relation of the acoustic mode 
frequency of an fee V-planes ultrathin film, described by 
Hamiltonian (1), is found to be, in the low wavevector 
limit 



n„(k||) « 7v^i(k||) Mh) - [ImB(k||)] 2 (Al) 
where 

yli(k||) = |ff [sin0sin0# cos(<,s — ipn) + cos9cos9h] 
+ H K cn cos 2 9 — H Tf± cos 4 9 

2 4 

+ -rH^w sin 2 0[3 cos 2 9 - cos 4^(3 + sin 2 0)1 ) 

+ [—H dip sm 2 ((p - ip k )j ■ (k\\a) 

+ H ex -(k\\af (A2a) 

A 2 (k||) = ao — d\ ■ (k\\a) + a 2 ■ (k\\a) 2 

= Ih [sin sin 0# cos^ — ipn ) + cos cos 9h] 
+ H K *t cos(20) - H K ± cos 2 0(1 - 4 sin 2 0) 

- -#„,! sin 2 0(1 -4cos 2 0)(3 + cos4v?)| 

4 4 J 

- [y H dip [sm 2 9 -cos 2 9 cos 2 (ip-<p k )}} 
■ ( k \\ a ) 

+ H ex -(k {l a) 2 (A2b) 

ImB(k||) = | — "^H K \i sin 2 0cos0sin4y>| 
rjV -I 

- <^—H d i P cos 9 sin[2(ip - cpk)]j ■ (k\\a) (A2c) 

The expression for fio(kii) is approximated in the sense 
that it is the generalizatirm of the spin-wave dispersion 
relation of the monolayeiOEjO to the case of an ultra- 
thin film with N planes£3 for the latter case, it repre- 
sents the energy of the acoustic mode of the film, which 
is measured in a BLS experiment. We notice that the 
effect of the finite number of planes manifests itself as 
a prefactor N in the terms containing H d i p and linear 
in the wavevector ku. In the previous equations, k|| is 
the in-plane two-dimensional wavevector forming an an- 
gle tpk with the in-plane x = [100] crystallographic axis; 
a = ao/v2 is the constant of the square lattice on the 
(001) surface expressed in terms of a , the constant of 
the bulk fee lattice; and ip denote the polar coordi- 
nates of the magnetization taking z (the normal to the 
film plane) as polar axis, while 9h and <pjj are the po- 
lar coordinates of the applied magnetic field. Dipolar 
sums were evaluated in the limit of small kii using an 



4 



Ewald-type summation method. 113 The dipolar field 
is Hdi P = 47rMoCi, where Mq = Ag^iBS/a^ is the mag- 
netization of the bulk fee lattiea-and c\ = f/V% is a 
thickness dependent coefficicnt£3e3 The other fields are 
defined in the text after Eqs. (|2a,2b). 

In the general case, the ground state configuration 
(9, ip) is obtained solving the system 



[H 
-\H 



-H„\\ (3 + cos Akp)] sin" 1 9 cos 9 

4 4 

H K ea \ sin 9 cos 9 



-H[cos9 sm9n cos(ip — ipjj) — sin 9 cos Ojj] — (A3a) 



H sih9h sbx((f — tp H \ — -JJ„y sin 3 6* sin 4^ = (A3b) 

2 K 4 

When the external magnetic field is applied in plane 
along the easy direction (9h = vr/2, tp H = 7r/4), Eq. 
( A3b) is satisfied by ip = 7r/4, while Eq. (A3a) becomes 



cos9[(H K i 



-(H K x-H 



off sin i 



-H] = 



(A4) 



Hence, the zero-field canting angle reported in Eq. (j^) 
is obtained, provided that H K x > H K e«. Moreover we 
observe that, for ip = ipn = tt/4, ImB(k||) reduces to a 
term linear in Hdip and in the wa veve ctor: thus its square 
can safely be neglected in Eq. (Al), so that Eq. (||) is 
obtained. Finally, it is worth noticing that in the case 
of zero in- plane anisotropy, iJ„n = 0, one recovers for 

^4 

Ai(k||) and A2(k||)_the expressions previously found by 
Erickson and Mills .0 
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It is worthwhile noticing that in Co/Cu multilayers a direct 
visualization of the domaki structure using magnetic force 
microscopy was possible.Ej 



FIG. 2. Brillouin light scattering data for the spin wave 
frequency fio of a t = 60 A-thick Ni film versus the intensity 
of a magnetic field H applied along the in-plane easy axis 
[110], on increasing H (solid squares) as well as on decreasing 
it (solid circles). The solid line is the spin wave frequency, 
calculated from Eq.(^) assuming a uniform in-plane magne- 
tization. In the inset we compare the Qo(H) data measured 
for H applied along the easy axis [110] (solid symbols) and 
for H along the hard axis [100] (open symbols). 



FIG. 1. Brillouin light scattering spectra for several val- 
ues of the external field H, applied in-plane along the [100] 
direction. The Ni film was 60 A thick, and data were ob- 
tained using the backscattering configuration; see text for the 
experimental details. 



FIG. 3. Magneto-optic vector magnetometry data in a 
canted field geometry, showing the field dependence of: a) 
the polar (M z ) component of the magnetization, b) the lon- 
gitudinal (M x ) component, c) the magnetization modulus 
M = V M'i + M'i . All quantities are normalized to the sat- 
uration value Msat- The field was applied at different angles 
(6 H = 0°, 5°, 11°, 20°) with the film normal, z. 
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